Astron. Nachr. / AN 327, No. 9, 844 -18491(2006) / DPI 10.1002/asna.200610662 



Nonlinear simulations of magnetic Taylor- Couette flow with current- 
free helical magnetic fields 

J. Szklarski* and G. Riidiger 

Astrophysikalisches Institut Potsdam, An der Stern warte 16, D- 14482 Potsdam, Germany 

Received July 2, 2006, accepted later 
Published online more later 

Key words methods: numerical - magnetic fields - magnetohydrodynamics (MHD) 

The magnetorotational instability (MRI) in cylindrical Taylor-Couette flow with external helical magnetic field is simulated 
for infinite and finite aspect ratios. We solve the MHD equations in their small Prandtl number limit and confirm with time- 
dependent nonlinear simulations that the additional toroidal component of the magnetic field reduces the critical Reynolds 
number from O(10 6 ) (axial field only) to O(10 3 ) for liquid metals with their small magnetic Prandtl number. Computing 
the saturated state we obtain velocity amplitudes which help designing proper experimental setups. Experiments with 
liquid gallium require axial field ~ 50 Gauss and axial current ~ 4 kA for the toroidal field. It is sufficient that the vertical 
velocity u z of the flow can be measured with a precision of 0. lmm/s. 

We also show that the endplates enclosing the cylinders do not destroy the traveling wave instability which can be observed 
as presented in earlier studies. For TC containers without and with endplates the angular momentum transport of the MRI 
instability is shown as to be outwards. 

© 2006 W1LEY-VCH Verlag GmbH & Co. KGaA, Weinheim 

1 Introduction 

The magnetorotational instability (MRI) has been formu- 
lated long time ago (Velikhov 1959) but recently it become 
of particular interest for astrophysics due to recognition 
of its importance as the source of turbulence in accretion 
disks with Keplerian rotation (Balbus & Hawley 1991). One 
of the most convenient models for MRI is magnetohydro- 
dynamical cylindrical Taylor-Couette flow, i.e. motion of 
liquid metal between two concentric differentially rotating 
cylinders under imposed external magnetic field along the 
axis (Riidiger & Zhang 2001; Ji, Goodman & Kageyama 
2001; Willis & Barenghi 2002). Under the influence of a 
purely axial field MRI operates for magnetic Reynolds num- 
ber Rm = Oi n i?i n (i? out — Rin)M of order O(10) and has 
never been observed in an experiment due to very small 
Prandtl number for laboratory liquids (see Table ^ which 
results in need of high critical Reynolds numbers of order 
10 6 . Recently it has been shown by Hollerbach & Riidiger 
(2005) and Riidiger et al. (2005) that a current-free external 
toroidal magnetic field in addition to the usual axial field can 
reduce the critical Reynolds number to O(10 3 ) and there- 
fore it makes it much easier to design a MRI experiment 
due to slower rotation. 

The aim of the presented paper is to compute the ampli- 
tudes of velocities and the frequencies in a nonlinear satu- 
rated state under the influence of the endplates covering the 
cylinders. 

The fluid confined between the cylinders is assumed to 
be incompressible and dissipative with the kinematic vis- 
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Fig. 1 The geometry of the problem - two concentric 
cylinders with radii R lu — 4 cm, i? ou t = 8 cm rotating 
with f2j n , f2 ou t- B z and are the external magnetic fields. 
Bcj, is due to an axial current inside the inner cylinder. In 
the PROMISE experiment cylinders have finite length of 
H = 40 cm (r = 10) and are covered with endplates. The 
liquid is a mixture of gallium, indium and tin (Stefani et al. 
2006). 



cosity v and the magnetic diffusivity rj. Derived from the 
conservation of angular momentum the rotation law Q(R) 
in the fluid is 



_b_ 

m 



(i) 
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Table 1 Physical properties of liquid metals which 
are suitable for the MRI experiment. A star denotes the 
PROMISE mixture. 





P [g/cm 3 ] 


v [cm^/s] 


r) [cm'Vs] 


Pm 


sodium 


0.9 


7.1 ■ 10" 3 


0.8 • 10 3 


0.9 • 10" 5 


gallium 


6.0 


3.2 ■ 10~ 3 


2.1 ■ 10 3 


1.5 ■ 10~ 6 


gallium* 


6.4 


3.4 ■ 10~ 3 


2.4 ■ 10 3 


1.4- 10~ 6 



with 



" — '., ^2 in m J 



(2) 



1 — ?7 Z 1 — i] 

where jl = f2 out /f2i n . After the Rayleigh stability criterion, 
i.e. dn(R 2 Vl) 2 > 0, the flow is hydrodynamically stable for 
ji > fj 2 , i.e. ji > 0.25 for our geometry. In order to isolate 
MRI we are mainly interested in hydrodynamically stable 
regimes so that ji > fj 2 should be fulfilled. 

2 Equations and numerical model 

The dimensionless incompressible MHD equations are 

Ha 2 

d t u + (u ■ V)it = -Vp + V 2 it H (rot B) x B, 

Pm 

d t B = —X7 2 B + rot(ux B), (3) 
Pm 

with div u = div B = and with the Hartmann number 



Ha = B, 



-Rin(-Rout — -Rin) 



(4) 



The magnetic Prandtl number is Pm = v/r). D = i? out — 
Rj n is used as the unit of length, v/D as the unit of velocity 
and B as the unit of magnetic fields. The Reynolds number 
Re is defined as Re = Vt m R ln D/v. 

For laboratory liquids like gallium, the magnetic Prandtl 
number is very small (see Tabled so the magnetic diffusion 
time is much shorter than other time scales and fluctuations 
b of the field Bq + b are also much smaller than the applied 
field Bo- As a result the b adjusts instantaneously to the 
velocity u and the quasi-static approximation can be used 
(Roberts 1967; Zikanov & Thess 1998; Youd & Barenghi 
2006). 

With B — > B + Pmb and Pm — > one finds the 
equation system 

d t u + (u ■ V)u = -Vp + V 2 tt + Ha 2 (V x b) x B , 



V 2 b = -rot(tt x B ) 

and div u = div 6 = 0. The external field is 



B = (0,B^B ) = B 



R 



1 



(6) 



where the parameter (3 denotes the ratio of toroidal field 
at the inner cylinder to the constant axial field Bo . Except 
along the axis the magnetic field is current-free. 

We simulate the axisymmetric 2D flow in cylindrical 
coordinates (R, </>, z) with the numerical code of A. Youd 



(see Youd & Barenghi 2006 for details). The code has been 
modified in order to handle periodic boundary conditions 
in a way suitable for our needs, the toroidal field was added 
and different boundary conditions on the endplates were ap- 
plied. 

Equations l|5} are solved with a finite-difference method 
in the R-z plane in a streamfunction-vorticity formulation, 
i.e. 

d tU<f> = (V 2 - ir 2 )u + [u x (rotu)] + Ha 2 <9 ZJ B 

dtuj^ = (V 2 - R- 2 )uj^ 

[rot (w x + Ha 2 [rot {j x B )]<j>, (7) 

with 'Poisson equations' for streamfunction, magnetic field 
and the current 



(- J R- 1 (9 zz + 9 rr )+ J R- 2 9 r )^ 
(V 2 B) 

/V7 2 -x a f \ , a 9r{Ru r ) 

(V 3)4, = -d zz (u R ) + d R - 



= ^0 

d z {u4, - u z f3B —^) + dRUR(3B —^ 
H K 



R 



(8) 



where the streamfunction ip is defined by u R = 
— (1/R)d z ip, u z = (l/R)d R 4> and vorticity u> = rot (it). 

No-slip boundary conditions for the velocity on the 
walls and the endplates are used, so that u r — at R = 
Rm, Rout and u z — at z = 0, H. Therefore ip = at 
R = Rj u , Rout and z = 0, H. Using ©i we obtain lo^ — 
-(1/R)d rr ip at R = i? in ,i?out and ui^ = -(l/R)d zz %p at 
z = Q,H. Boundary conditions for u^, at the walls are de- 
termined by Re and ji and at the endplates by their rotation 
properties - fixed or rotating as rigid body with f2 out or f2j n . 

The magnetic boundary conditions depend on the elec- 
trical properties of the walls. For perfectly conducting cylin- 
ders we assume infinite conductivity so that with J = 
<t(E + it x b) and E$ — E z — for the walls we get 

R-% + dub^ = 0, j = O at R = R in ,R out , (9) 

In the present paper only insulating endplates are considered 
and the BCs for magnetic field are obtained by applying 
the so called pseudo-vacuum approximation for which b r = 

= d z b z 

= 9 2 J0 

at z = 0, H. We use grid resolution Nr. x N z = 8Q x 800 
(and 40 x 400 for comparison, the difference in obtained 
values is always less than couple of percent). The time-step 
is constant with dt = 10~ 5 or dt = 10~ 4 . 



0. Then 





(10) 



(5) 3 Infinite container 



First we have performed simulations for periodic verti- 
cal boundary conditions in order to compare the critical 
Reynolds numbers with those presented in Riidiger et al. 
(2005) as well as to compute for the first time amplitudes of 
velocities in the nonlinear saturated state. For infinite cylin- 
ders and jl = 0.27 the flow is always hydrodynamically 
stable, with external axial magnetic field it looses the sta- 
bility for large Reynolds number of order O(10 6 ) but the 
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Fig. 2 Snapshots of contour lines of the streamfunction: 
dashed lines represent counterclockwise rotation, solid lines 
clockwise rotation. Ha = 9.5, Rc = 900, fi = 0.27, = 4. 
The graph is not at scale. 
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Fig. 3 Frequency of the traveling wave as a function of 
rotation ratio fi for infinite cylinders, Ha = 9.5, Re = 
900, j3 = 4. Here and below the units are scaled according 
to physical properties of gallium. 
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Fig. 4 Amplitudes of velocity as function of Reynolds 
number, u'^ is the deviation from the standard Couette solu- 
tion. The amplitudes increase from the critical point Rc c = 
842 and continue to increase with Re. Infinite cylinders, 
Ha = 9.5, fi = 0.27, /? = 4. 



3.1 Frequencies 

Figure[2]shows snapshots of the streamfunction for the satu- 
rated state for Reynolds number just above the critical value, 
the three pairs of drifting vortices can be seen. This drift 
and its velocity agrees with the previous results, it is slow, 
for gallium the period of the traveling wave is about 230 s 
for j3 = 4 which corresponds to Wdrift = 0.5 mm/s (the pe- 
riod of one rotation of the inner cylinder is about 35 s for 
Re = 900). 

From the experimental point of view it is also conve- 
nient to present traveling wave frequencies as function of fi 
(Fig. [3J. All the frequencies were obtained by taking the 
spectrum of the velocity component u z (t) in the middle 
of the gap and then choosing the dominant frequency. Ex- 
cept for very small ft the resulting frequencies are always 
smaller than the frequencies of both the cylinders. Beyond 
the Rayleigh limit (i.e. for fi > 0.25) the characteristic fre- 
quency is about 10% of the frequency of the inner cylinder. 
That means that beyond the timescale of rotation another 
timescale exists in magnetic astrophysical systems which 
exceeds the rotation period by (say) a factor of ten. Due 
to the small electric conductivity, however, our TC-system 
only possesses a very short diffusion time of order of 10 _2 s. 



additional toroidal current-free field can reduce this number 
to O(10 3 ) -for j3 = 3, 4 the critical Re crit is 1160 and 842, 
resp. 

In the simulations the length of the periodic cylinders 
was chosen to be three times the wavelength obtained with 
linear analysis. The agreement between the previous and 
presented results is rather good. 



3.2 Velocity amplitudes and torque 

The saturated nonlinear solution provides the values of ve- 
locity amplitudes which have not been known before. The 
velocity component which is easiest to measure in the labo- 
ratory is u z . 

Figure|4]shows the dependence of different components 
of velocity as function of Reynolds number Rc. The veloc- 
ities are zero up to critical point Rc cl -i t at which they start 
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, 5 Average maximum amplitudes of the velocity ver- 
fi for infinite cylinders, Re = 900, (3 — 4, Ha = 9.5. 
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Fig. 6 Torque at the walls for two different (3. Note that 
from the critical Rc cl -i t = 842 for (3 = 4 case the torque is 
larger than for the stable (3 = case, Ha = 9.5. 



to grow. We notice that the values of u can be further in- 
creased by increasing the rotation rates of the container and 
therefore it is easier to measure u z for higher Re. 

In Fig.[3]the maximum amplitudes of velocity for Re = 
900 and (3 = 4 are presented. For any stable solution all the 
amplitudes should be zero just beyond the Rayleigh line, 
i.e. on the right of /t = 0.25. We see that with additional 
toroidal field this is not the case - the velocities for fi = 0.27 
are of order 0.1 mm/s and drop here to zero for fi = 0.28 
for these parameters. 

For MRI one expects the angular momentum to be trans- 
ported outwards and this should result in an increase of the 
torque. We are mainly interested in torques at the walls at 
which the magnetic and advective torque vanish so the total 
torque is equal to the viscous one, which is defined as 



Although the change of the torque is rather small and would 
be difficult to measure in the laboratory, Fig. [6] shows that 
increase of the torque can be observed in our computations 
above the critical Reynolds number. Indeed the angular mo- 
mentum is transported outwards. 

4 Finite container 

In any real experiment finite cylinders enclosed by end- 
plates are used. It has been pointed out that the endplates 
can have significant impact on the flow especially for short 
cylinders and for fast rotation (e.g. Hollerbach & Fournier 
2004; Kageyama et al. 2004) also some solutions for re- 
duction of the endplates effect have been proposed (Burin 
et al. 2006). We are mainly interested in knowing how the 
endplates and their mechanical and magnetic properties in- 
fluence the velocity amplitudes and frequencies for slow ro- 
tating finite cylinders with aspect ratio of (say) T = 10 and 
a current-free helical external field geometry. 

4.1 Fkman layer 

The existence of the viscous endplates results in developing 
so-called Ekman layers in which the velocity differs from 
any toroidal rotation as centrifugal force is weakened close 
to the vertical boundaries. This effect drives a global merid- 
ional circulation and two Ekman vortices appear. The veloc- 
ity in this region depends on the Reynolds number and the 
velocity of the endplates themselves. As known for nonro- 
tating endplates (or rotating as a solid body with 5! out ) there 
is a radial inflow close to the boundaries and for solid-body 
rotation with Q- m there is a radial outflow. 

In the Rayleigh-stable regime the flow is hydrody- 
namically stable for both the infinite and finite contain- 
ers, whereas for the latter two Ekman vortices are always 
present. We confirmed it with the numerical method that it 
is used here. 

A simple estimate of the thickness of the Ekman layer 
(for endplates rigidly rotating with £l ut, see e.g. Ji et al. 
2004) gives 



0.13 cm 



(12) 



T(R) 



H 



R 3 pvd R 



R 



dz. 



for the (3 = 4 critical case where fl = ^/fli n il out . This is 
a larger value in comparison to other experiments dealing 
with faster rotation (Re « 10 6 ) in which case the possibil- 
ity of developing hydrodynamical turbulence in the Ekman 
layer arises. 

The typical velocity in the Ekman flow can be approxi- 
mated by us = \J n^l which for Re = 900 with 0.02 cm/s 
is of the same order for velocities for the unstable case pre- 
sented for infinite cylinders (see Fig.|5Jl. For Re = 900 the 
maximum value of u z (in the corners) is above 0.1 cm s _1 
which is ten times more than the values computed for infi- 
nite cylinders. Moreover, for the same parameters the Ek- 
^ ' man circulation time H/2V v£l is of order 10 3 s which is 
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Fig. 7 The same as in Fig. [5] but for finite cylinders. 
Re = 900, (3 = 4, Ha = 9.5, insulating endplates (pseudo- 
vacuum boundary conditions) both rotating with the outer 
cylinder. 
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Fig. 8 Frequencies as function of rotation ratio fi for finite 
cylinders, a) (3 — 6, Re = 1500, Ha = 9.52, insulating 
endplates, bottom: fixed, upper: rotating with il ou t- b) (3 = 
4, Re = 900, Ha = 9.5, insulating endplates rotating with 

^out ■ 

comparable to typical MRI growth time. We must conclude 
that endplates can significantly alter the flow. 

In the computations it turns out that in the presence of 
the external axial field, the Ekman flow induces significant 
z-gradients in b<f, which can result in a magnetic instabil- 
ity even for small Reynolds number. The critical numbers 
depend on magnetic properties of the endplates, their rota- 
tion and Hartmann number 1 . Since at this moment our aim 
is to investigate how the traveling-wave effects of the MRI 
with helical (current-free) field structure can be observed 
experimentally, this instability, although highly interesting, 
is beyond the scope of the present paper. 

1 e.g. Re cr ; t = f40 for Ha = 3.5, perfectly conducting, fixed end- 
plates 
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P 

Fig. 9 Frequencies of the traveling waves in the middle of 
the gap for finite cylinders as function of (3. Re = 900, fi = 
0.27, Ha = 9.5, insulating endplates rotating with O out . 

4.2 Amplitudes and frequencies 

The velocity amplitudes for finite cylinders just above the 
critical Reynolds number are shown in Fig. [7] We notice 
that due to the Ekman circulation those velocities do not 
drop to zero even for (i = 0.28, which was true for infinite 
cylinders. Again the external toroidal field maintains trav- 
eling waves which can easily be observed in the laboratory. 
Analog to Fig. [5] Fig. [8] shows the frequencies of the ve- 
locity field in the center of the gap. Unlike for the infinite 
container, here the position of measurement does make dif- 
ference since in the vicinity of the endplates the traveling 
wave is suppressed by the Ekman layer. From Fig.[9]we can 
clearly see that non-zero j3 introduces vertical movement 
with very similar frequencies as for the infinite case (com- 
pare with /ft = 0.27 in Fig.[3j. 

Snapshots of streamfunction are shown in Fig. ^3 We 
see that the flow is complicated and deformed when com- 
pared with infinite-container solutions. The z-asymmetry 
close to the endplates is due to different boundary condi- 
tions: the bottom endplate is at rest, the upper moves with 
the outer cylinder 2 . We conclude that even for nonsymmet- 
ric rotation of the endplates the traveling wave can still be 
observed. 

Similarly like for the infinite case one can also observe 
the increase of torque in comparison with (3=0 case 
(Fig.Qj). 

5 Conclusions 

We have shown that, as predicted, additional toroidal fields 
reduce the critical Reynolds number for the MHD Taylor- 
Couette flows for laboratory liquids with Pm as 10~ 5 — 
10~ 6 . The nonlinear, axisymmetric simulations in the quasi- 
static approximation yield very similar characteristic values 

2 This is the case for the current experimental setup of PROMISE 
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Fig. 10 Snapshots of contour lines of the streamfunction 
analog to Fig.[2]but for finite cylinders. Rc = 1500, (3 — 
6, Ha = 9.5, insulating endplates, bottom: fixed, upper: ro- 
tating with f2 out . 
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Fig. 11 Torque at the walls (inner and outer) for finite 
cylinders for different values of j3. Re = 900, Ha = 9.5. 



compared to those obtained with linear theory for infinite 
cylinders. 

The velocity amplitudes which were not known until 
now, are small but still it is not very difficult to measure 
them experimentally. Similarly the velocity of the drift can 
easily be detected by looking at changes in the velocity field. 
We shall also mention that increasing Reynolds numbers re- 
sult in linear (at least to some point) increase of velocity 
amplitudes, therefore it is easier to detect them for higher 
Re, however for too high rotation rates the assumption of 
axisymmetry may not be valid. 



For finite aspect ratio calculations, the velocities in- 
duced by Ekman pumping are of similar order as those re- 
sulting from the instability. Therefore one should not ne- 
glect the influence of the endplates in such experiments. 
Nevertheless the influence of the toroidal field and the trav- 
eling wave can easily be detected in all cases. 
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